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ABSTRACT 

"ih : 

^ ! We have observed a high-mass protobinary system in the hot core W3(H20) 

I with the BIMA Array. Our continuum maps at wavelengths of 1.4 mm and 

^ ! 2.8 mm both achieve sub-arcsecond angular resolutions and show a double-peaked 

^ ! morphology. The angular separation of the two sources is I'.'IQ corresponding to 

^ ! 2.43 X 10^ AU at the source distance of 2.04 kpc. The flux densities of the two 

^ ! sources at 1.4 mm and 2.8 mm have a spectral index of 3, translating to an 

^ i opacity law of cc u. The small spectral indices suggest that grain growth has 

^ i begun in the hot core. We have also observed 5 K components of the methyl 

Q\ ', cyanide (CH3CN) J = 12 — > 11 transitions. A radial velocity difference of 

^ i 2.81 ±0.10 kms~^ is found towards the two continuum peaks. Interpreting these 

T— I I two sources as binary components in orbit about one another, we find a minimum 

^ [ mass of 22 Mq for the system. Radiative transfer models are constructed to 

^ I explain both the continuum and methyl cyanide line observations of each source. 

0-(\ Power-law distributions of both density and temperature are derived. Density 

O I distributions close to the free-fall value, r~^-^, are found for both components, 

I suggesting continuing accretion. The derived luminosities suggest the two sources 

I have equivalent zero-age main sequence (ZAMS) spectral type BO. 5 - BO. The 

I nebular masses derived from the continuum observations are about 5 M0 for 

^ I source A and 4 M© for source C. A velocity gradient previously detected may be 

I explained by unresolved binary rotation with a small velocity difference. 

Subject headings: circumstellar matter — stars: individual (W3(0H), W3(H20)) 
— stars: pre-main sequence — stars: early-type 
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1. INTRODUCTION 

Massive star formation presents a challenge in both observational and theoretical stud- 
ies. Unlike their low-mass counterparts, massive young stellar objects (YSOs) are typically 
distant {d > 0.5 kpc), obscured {Ay > 100), and often observed in clustered environments. 
With the development of millimeter interferometers, it has become possible to resolve indi- 
vidual, deeply embedded, massive YSOs. In the past two decades, a few warm {T > 100 K), 
dense {n > 10^ cm~^), molecular clumps of high luminosities (L ~ 10^ — 10^ L©) have 
been discovered in the vicinity of ultracompact (UC) H ll regions (Turner & Welch 1984; 
Kurtz et al. 2000). These molecular clumps exhibit chemistry that is characterized by high 
abundances of fully hydrogenated species (van Dishoeck & Blake 1998), including a variety 
of organic molecules, such as methyl cyanide (CH3CN), ethanol (C2H5OH), etc. In a few 
cases, internal heating of the clumps is showed directly by radial temperature gradients with 
ammonia (NII3) emission (Cesaroni et al. 1998). These molecular clumps were later classi- 
fied as hot cores. The lack of detectable Stromgren spheres further places hot cores in an 
evolutionary phase earher than UC Hll regions (Wilner et al. 2001). Maser emission, such 
as water and methanol masers, is often associated with hot cores and suggestive of violent 
collision events (Tarter & Welch 1986). One most interesting phenomenon is an apparent 
velocity gradient that has been detected in a few hot cores, for example, the Orion hot core, 
W3(H20), IRAS 20126+4104, NGC 7538 S, G29.96-0.02, G31.41+0.31, and IRAS 07427- 
2400 (Plambeck, Wright, & Carlstrom 1990; Wyrowski et al. 1997; Cesaroni et al. 1999; 
Sandell, Wright, & Forster 2003; Olmi et al. 2003; Cesaroni et al. 1994; Kumar et al. 2003). 
Although the emission is spatially unresolved in each velocity channel, the velocity gradient 
has been long proposed as evidence for the Keplerian rotation of a circumstellar disk. 

At a distance of 2.04 kpc (Hachisuka et al. 2006), the neighborhood of the W3(0H) 
ultracompact (UC) H ll region is a nearby, well-studied star-forming region in the course of 
developing an OB stellar group. The most recent star-forming activity is taking place in a hot 
molecular core, W3(H20), which appears as a warm, dense molecular clump 6" to the east of 
W3(0H). This hot core was first discovered by Turner & Welch (1984) through an aperture 
synthesis study of HCN (1 — > 0) emission. The broad line width observed in HCN indicated 
the presence of an obscured young stellar object, possibly a protostar, with an estimated 
luminosity of 10^ L0. Proper motions of water masers in W3(Il20) can be explained as 
part of a bipolar outflow (Alcolea et al. 1993). Wilner, Welch, & Forster (1995) detected 
unresolved dust emission at 87.7 GHz with a flux density of 43 mJy, implying a mass of 
10-20 Mq and a luminosity of 10^-10^ L©. At centimeter wavelengths, observations with the 
Very Large Array (VLA) revealed a synchrotron jet near the expansion center of the water 
masers (Reid et al. 1995; Wilner, Reid, & Menten 1999). A secondary continuum source with 
a rising spectral index was also detected, possibly featuring an additional young stellar object 
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(Wilner et al. 1999). Subsequent line and continuum observations at 220 GHz (Wyrowski 
et al. 1999) resolved the hot core into three emission components, with two components 
(A and C) roughly coincident with two temperature peaks of about 200 K, derived from 
the HNCO Ka = 0, 2, and 3 lines. In addition, a linear velocity gradient spanning about 
10 kms~^ along the East-West direction was found within a region of 1" across the hot core 
(Wyrowski et al. 1997). This gradient was found by plotting the positions of maximum 
emission derived from channel maps in several lines. The total gradient spanned an angular 
range of about 1", smaller than the instrumental resolution. Thus, it was unclear what 
could be the source of the gradient. Wyrowski et al. (1997) proposed three interesting 
possibihties. The first suggestion was that the high excitation molecules traced the sides of 
the cavity carved out by the outflows associated with the water masers. The second idea was 
that the gradient represents a large disk seen edge-on. The third idea was that there might 
be a double source with different radial velocities with the components associated with the 
two centimeter continuum sources. 

The present study is an extension of the Plateau de Bure work (Wyrowski et al. 1997, 
1999) enabled by the Berkeley-Illinois-Maryland Association (BIMA) Millimeter Array with 
its longer baselines and more extensive uv coverage. Although the BIMA Array has less 
collecting area, the W3(H20) source is sufficiently bright to be imaged. We have mapped the 
continuum emission at wavelengths of 1.4 mm and 2.8 mm at very high angular resolution 
in search of the large disk or whatever else might be the source of the velocity gradient. 
The angular resolution is 0'/26 at 1.4 mm and 0'/40 at 2.8 mm, allowing us to study the 
structure and dust properties of the components within the hot core in detail. We have also 
observed the K = 2,3, 4, 5, and 6 components of the CH3CN J = 12 — 11 transitions with 
1" resolution to derive the temperature at each continuum peak. The kinematics are also 
obtained with high precision by simultaneously fitting Gaussians to each of the 5 observed 
K components. 



2. OBSERVATIONS AND DATA REDUCTION 

The continuum emission was mapped at both 1.4 mm and 2.8 mm in the A and B config- 
urations of the BIMA Array The pointing center was (a,(5)(J2000) = (2^'27"'3!87, +6r52'24'.'6) 
The calibrators were 0359+509 at 1.4 mm and 0228+673 at 2.8 mm. The flux densities of the 
calibrators were determined by comparison with the Be star MWC 349, whose flux densities 
were assumed to be 1.2 Jy and 1.9 Jy at 2.8 mm and 1.4 mm, respectively (Dreher & Welch 
1983; Tafoya, Gomez, & Rodriguez 2004). The accuracy of the flux scale was estimated with 
the variance in the flux measurements of 3C84, and should be good within 15%. 
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The 1.4 mm data set taken in the A configuration provided highest spatial frequencies, 
which resolved the UC Hll region W3(0H) into a shell-like structure. The three continuum 
data sets at 1.4 mm were observed in fast-switching mode, consisting of a 2 minute cycle 
on 0359-1-509, 0228-1-673, and W3(0H) to trace the rapid atmospheric phase fluctuation. 
Additionally, a 1.5 minute integration on 0359-1-509 was repeated every 20 minutes to help 
gain amplitude solutions. The reference source, 0228-1-673, was included to monitor the 
imaging quality. The typical system temperature was 400 ~ 600 K, and the typical phase 
noise was about 0rms — 35°. The 2.8 mm continuum data sets were observed in a regular 
fashion, which repeats the calibrator, 0228-1-673, every 5 minutes. 

The calibration and imaging were performed with the MIRIAD software package. The 
two sidebands were calibrated independently. Every data set was first calibrated by deriving 
phase-only gain solutions from the calibrator on a timescale of 2 or 5 minutes, followed by 
solving complex gains on a longer timescale of 15 to 20 minutes. The combined complex 
gains were transfered to the source data, which were further improved with one iteration 
of phase-only self calibration, with the exception of the 1.4 mm A-configuration data that 
showed no improvement after self calibration. The self-calibration intervals were about 1 
or 3 minutes, depending on the wavelengths and configurations. The self calibration was 
performed on W3(0H) for the A-configuration data and on both W3(0H) and W3(H20) for 
the B-configuration data. When multiple data sets were available, a combined map was made 
to be the model for self calibration. The astrometery at 1.4 mm after the self-cahbration 
was checked with reference quasar, 0228-1-673: first apphed the phase-only gain solutions 
to 0228-1-673, generated an image, and then performed Gaussian fit to find the emission 
ccntroid. The positional errors were less than 0'.'014 and 0'.'05 for the A- and B-configuration 
data, respectively. Because the observations in the A configuration resolved W3(H20), the 
images were formed by combining the visibilities from both the A and B configurations. Both 
sidebands were used in a mult i- frequency synthesis to generate line-free continuum maps. 
Each map was then CLEANed in the usual way and was restored with a circular Gaussian 
beam of the same sohd angle as its synthesized beam. Details of the continuum observations 
are listed in Table 1. 

We have also observed the = 2, 3, 4, 5, and 6 components of the CH3CN J = 12 — > 11 
transitions in the B configuration during the spring season of 2003. In order to improve 
the single-to-noise ratio (SNR) on the calibrator, 0359-1-509, for each 5 minutes interval, 
the correlators were adjusted to include two wide-band windows, which provided 400 MHz 
continuum bandwidth. However, this correlator setup allowed simultaneous integration of 
only two high spectral resolution windows. The K — 2,3, and 4 components were observed 
on February 8 while the K = 5 and 6 components were observed on February 9 under nearly 
identical weather conditions. The data were reduced with a procedure similar to that of the 
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Table 1. Parameters for the Continuum Observations 



Parameter 


1.4 mm 


2.8 mm 


Synthesized frequency 


221.2 GHz 


108.6 GHz 


Dates (configurations) 


2000 Dec 26 (A) 


2003 Jan 8 (A) 




2001 Jan 31 (B) 


2003 Jan 16 (A) 




2001 Feb 8 (B) 


2003 Feb 2 (B) 


Cahbrator 


0359+509 


0228+673 


uv-range (kA) 


53 - 1000 (A) 


16 - 460 (A) 




10 - 170 (B) 


4 - 85 (B) 


Largest visible scale 


2" (A) & 10" (B) 


6" (A) & 25" (B) 


W3(0H) flux density 


3.41 Jy 


3.61 Jy 


Primary beam 


52" 


106" 


Weight Scheme 


uniform 


uniform 


Synthesized beam 


0'/31 X 0'/22 (77°) 


0'/46 X 0'/35 (72°) 


Restored beam 


0"26 


0"40 


Bandwidth (each sideband) 


800 MHz 


400 MHz 


Rms noise (mJy/Beam) 


3.0 


1.7 
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continuum data. The velocity resolution of the spectral hne observations is 0.53 kms ^. The 
line observation parameters are listed in Table 2. 

3. RESULTS 

3.1. Continuum emission 

Fig. 1 shows the millimeter continuum emission (contours) for both 1.4 mm and 2.8 mm 
in the central part of the W3(0H) region overlaid on a 3.6 cm continuum map (Wilner et al. 
1999, greyscale). The brightest source is the UC Hll region W3(0H), whose champagne flow 
to the northeast is visible at 3.6 cm. The hot core W3(H20), 6" east to W3(0H), is detected 
at both 1.4 mm and 2.8 mm with flux densities of 1.38 Jy and 0.22 Jy, respectively. The 
dust emission is much stronger at 1.4 mm. In addition, a thin layer of dust may still wrap 
around the UC H ll region and contributes to the weak emission feature along the edge of 
W3(0H) in the 1.4 mm map. We do not detect other sources of continuum emission within 
the held of view. 

For the components in the hot core, we adopt the notation of Wyrowski et al. (1999), 
who identifled three peaks: A, B, and C. The principal features are A and C in our map. 
This is the flrst time that the separate source C has been imaged at 2.8 mm. The B feature 
is roughly the northeast extension of peak C in our 1.4 mm map. The 1.4 mm map of Fig. 1 
agrees best with the 3.6 cm map (Wilner et al. 1999). The peak positions are at offset A 
= (5'.'94, O'.'OS) and C = (4'.'75, O'.'IO) relative to the pointing center. This gives a projected 
separation of r.'19, equivalent to 2.43 x 10^ AU. 

The hot core features of Fig. 1 are best interpreted as two separate sources, possibly a 
binary system. Source B may be a separate companion to source C. In any case, a large disk 
seems very unlikely. There is very little emission bridging sources A and C, not even enough 
to flt an edge-on doughnut. Thus, we adopt the third idea of Wyrowski et al. (1997) that 
the structure of the hot core corresponds to two principal components, and we develop our 
models for the region with the idea that we are observing a protobinary system. In Sect. 5 
and 6, we show how this two component model nicely explains the large velocity gradient 
discovered by Wyrowski et al. (1997). Indeed the small deviations from the straight hne 
that are evident in our data are difficult to explain with any other model. 
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Fig. 1. — BIMA high-resolution hne-free continuum maps at 1.4 mm (contours, top panel) 
and 2.8 mm (contours, bottom panel) overlaid on a 3.6 cm continuum map (greyscale, Wilner 
et al. (1999)). The brightest source is the ultracompact Hll region, W3(0H). The hot core 
W3(H20), located 6" east to W3(0H), shows a double-peaked morphology at both 1.4 mm 
and 2.8 mm with peak positions coinciding well with those in the 3.6 cm map. The three 
crosses indicate the positions of the three continuum peaks. A, B, and C, from East to 
West, observed by Wyrowski et al. (1999). The small bar at lower-left corner indicates 
10^ AU. The beam size is 0'.'26 for the 1.4 mm map and 0'.'40 for the 2.8 mm map. Contour 
levels of the 1.4 mm map correspond to -12, -6, 6, 12 to 96 by 6 {2a), and 96 to 222 by 12 
(4o") mJy/Beam. The white contours in W3(0H) start at 30 mJy (lOcr), which is described 
by light grey contours in W3(H20). Contours of the 2.8 mm map correspond to -6.8, -3.4, 
3.4 to 23.8 by 3.4 {2a), and 23.8 to 533.8 by 51 (30cr) mJy/Beam. The white contours start 
at 23.8 mJy (14cr), which is the highest contour level in W3(H20). 
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3.2. Methyl cyanide line emission 

The K = 2,3,4,5, and 6 components of the CH3CN J = 12 ^ 11 transition were ob- 
served with a 1" beam. The parameters of the observed K components are hsted in Table 3. 
We did not observe the K = and 1 components because their frequency separation of 
4.25 MHz is not sufficiently large to avoid overlap with each other. The large hne width 
expected in star-forming regions (> 2 MHz in the case of W3(H20)) makes these two com- 
ponents largely unusable for most purposes. Fig. 2 shows the channel maps of the K — 3 
component. The UC Hll region, W3(0H), centered at (0,0), also shows methyl cyanide emis- 
sion between Vr = —49 to —45 kms~^. The two triangles indicate the positions of the two 
continuum peaks, A and C. Although each of the component A and C remains unresolved, 
the hot core is well resolved along the East-West direction. The centroids of the methyl 
cyanide emission show a position drift from East to West with increasing radial velocity. 
The angular resolution of 1" is smaller than the separation of l'/19 between source A and 
C. This allows us to obtain the spectrum for one source without much interference from the 
other source. The properties of methyl cyanide molecules are discussed in Appendix A. 



4. RADIATIVE TRANSFER MODELS 

In order to derive the density and temperature profiles, spherically symmetric models 
are developed for self-consistency to explain the continuum and methyl cyanide line emis- 
sions of each source A and C. The density and temperature are assumed to have power- law 
dependences on radius, 
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n = no ( — ) cm"^ (2) 
T = To ( - ) ' K, (3) 



where and Tq are the mass density of dust, number density of methyl cyanide molecules, 

and temperature at the reference radius rg = 500 AU, respectively. The gas and dust 
components are assumed to have the same temperature, Tg = = T, everywhere in our 
models. This assumption will be discussed in detail in Sect. 6.6. 

Model images which are convolved with the observed angular resolutions are used to 
compare with the maps. First, intensity images are made by integrating the intensity, 1^, 
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Table 2. Parameters for the CH3CN Line Observations 

Parameters 2003 Feb 08 2003 Feb 09 

Configuration B B 

Calibrator 0359+509 0359+509 

Observed K components 2,3,4 5,6 

uv range (kA) 14 - 175 14 - 175 

Largest visible scale 7" 7" 

W3(0H) flux density (Jy) 3.57 3.45 

Syntfiesized beam 1'.'2 x 0'.'9 (-62°) l'/2 x 0'/9 (-63°) 

Restored beam I'/O I'.'O 

Channel width 0.53 km s~^ 0.53 km s~^ 

Channel rms 2.9 K {K = 2, 3) 3.0 K (if = 5) 

2.6 K{K = A) 3.0 K (X = 6) 



Table 3. Observed K Components of CH3CN ( J = 12 - 11) Transitions 



K 


Frequency (GHz) 


giK ^ 


SjK ^ 


QikSik 


EjK (K) ^ 


Einstein A {s ^) 


2 


220.730266 


6 


140/12 


70.0 


97.4 


8.98 X 10"^ 


3 


220.709024 


12 


135/12 


135.0 


133.2 


8.66 X 10-^ 


4 


220.679297 


6 


128/12 


64.0 


183.1 


8.21 X 10-^ 


5 


220.641096 


6 


119/12 


59.5 


247.4 


7.63 X 10-^ 


6 


220.594438 


12 


108/12 


108.0 


325.9 


6.92 X 10-^ 



^ Statistical weight. 
^ Line strength. 
Upper level energy. 
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50 50 50 50 

RA offset (arcsec) RA offset (arcsec) RA offset (arcsec) RA offset (arcsec) 



Fig. 2. — Channel maps of the K = 3 component of the CH3CN J = 12 — 11 transition. The 
two triangles indicate the two continuum peaks, A and C. In the hot core, the centroid of the 
emission clearly shows a position drift from East to West with increasing radial velocity. The 
UC Hll region, W3(0H), is centered at (0,0) with the CH3CN emission between Vr = —48 
to —43 kms~^. Contour levels start from -5.7, 5.7, 10.4 to 57 by 5.7 K (2cr). 
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along the line of sight (LOS) with the radiative transfer equation 




(4) 



where is the source function and r. 



J aydz is the optical depth, which is an integral of 



the absorption coefficient a^, along the LOS. The Boltzmann level population is assumed for 
the methyl cyanide lines (see the discussion in Sect. 6.2). Hence the source function is the 
Planck function, — IkTv^ jc?^ for both the continuum and methyl cyanide line models. 
The intensity images are then convolved with the restored Gaussian beam to generate model 
images that can better simulate the observations. In order to select the model that best 
describes the observation results, the chi-square, x^; is calculated and optimized with the 
Levenberg-Marquardt method. The goodness of the fits is evaluated by the reduced chi- 
square, = x^//) where / is the degrees of freedom in the fits. Since the expectation value 
of the chi-square is equal to the degrees of freedom, x^ — /, it is easier to compare the results 
with the reduced chi-square, whose expectation value is unity, — ^, ior a. moderately good 



Since our methyl cyanide observations did not resolve each individual component, we 
have chosen two power-law indices, q = 0.4 and 0.75, corresponding to the optically thin 
and thick regimes to the infrared (IR) photons, to investigate the density profiles and other 
physical quantities. We will explain how we determine the range of q in Sect. 5.2. The dust 
distribution is assumed to be a thick shell-like structure with an inner radius i?cav, which may 
correspond to the dust destruction radius, and an outer radius Rd, beyond which no dust is 
present. Four parameters are fitted for the continuum models: po, p, Rd, and i?cav) while q 
and To are held fixed. The methyl cyanide line models are also fitted with four parameters: 
uq, To, the Gaussian line width, b, and the radial extent of methyl cyanide, -RcHsCN- Again, 
two parameters, q and -Rcav, of line models are held unchanged during optimization. The 
line profile, (p^, is assumed to be a Gaussian function with a half- width, b, to 1/e: 



Lastly, the continuum and line models are iterated until self-consistency is obtained for the 
three parameters, Tq, p, and i?cav that are used in both models. 

Considering that molecules can not survive in the stellar ultraviolet (UV) radiation 
without shielding by dust grains, the sizes of the central cavities, Rcav, are assumed to be 
the same for the continuum and line models. Nevertheless, the derived temperatures are 
insensitive to the sizes of i?cav due to the large beam size of 1". The values of Tq do not 
fluctuate more than 5 K when i?rav increases from 20 to 300 AU. 



fit. 



1 ( "-"0 



(5) 
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5. CONTINUUM OBSERVATION RESULTS 

5.1. Grain growth in W3(H20) 

The continuum spectra of the two sources from 1.6 GHz to 221.2 GHz are shown in 
Fig. 3. The flux densities at 1.4 mm and 2.8 mm are found from this study by integrating 
the flux in a circular region of &!5 radius centered at each peak. The results are hsted in 
Table 4. The 15% systematic errors are included in the uncertainties. The spectral indices 
arc 3.0 ± 0.4 for source A and 2.9 ± 0.4 for source C. In summary, both sources show a rising 
spectrum between 1.4 mm and 2.8 mm with a spectral index of approximately 3.0. The 
spectral indices wc observed here arc smaller than, though consistent with, a spectral index 
of 3.6 in previous studies, which were based on a weak measurement at 3.4 mm (Wilner 
& Welch 1994). The improved sensitivity and angular resolution give us more reliable flux 
density measurements and the frequency dependence of the opacity can be determined for 
both sources. 

If the dust thermal emission is optically thin, containing a factor i/^ from the Planck 
function, a millimeter spectrum of F^, oc i/3.o±o.4 .^j^ ^^ve a mass opacity law of Ki, oc 
with P — 1.0 ± 0.4. On the other hand, a smaller spectral index could be a result of dust 
emission being optically thick at 1.4 mm. We will use our continuum models to show that 
the averaged optical depths are less than (r) ~ 0.09 in the cores (see Sect. 5.4 and Table 5). 
Hence, the small spectral indices are not the result of large optical depths. 

The frequency dependence of the dust opacity Ki, at long wavelength is an important 
parameter that can distinguish the particle types such as their shapes and mixtures. Dust 
grains in various environments may differ dramatically. The diffuse interstellar medium, 
consisting a mixture of small silicate and graphite spheroids, has a dust opacity spectral 
index of 2 at wavelengths longer than 100 fim (Draine & Lee 1984). In circumstellar disks, 
random aggregation of small dust grains may result in the formation of large particles with 
fractal dimensions. When grain growth happens, the emissivity at long wavelengths will be 



Table 4. Millimeter Continuum Flux Densities and Spectra 





Coordinate (J2000) 


FlA mm 


-^2.8 mm 




Source 


RA Dec 


(mJy) 


(mJy) 




A 


2^27°^4!709 +6r52'24'.'64 


474 ± 71 


55 ±8 


^3.0±0.4 


C 


2^27™4!541 +61°52'24^70 


351 ± 53 


44 ±7 


j^2.9±0.4 
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enhanced (relative to that at short wavelengths) and the spectral index of the opacity may 
decrease down to < 1 (Wright 1987; Miyake & Nakagawa 1993). A smaller spectral index 
at millimeter wavelengths has been proposed as evidence of grain growth in dense cores and 
in the circumstellar disks around low-mass pre-main-sequence stars (Beckwith et al. 1990; 
Chandler et al. 1995; Beuther, Schilke, & Wyrowski 2004). In a high density environment, 
collision and sticking can be much more effective to encourage grain growth. The observed 
dust opacity law of Ki, oc i/i o^o-^ suggests possible grain growth in the vicinities of the massive 
protobinary system in W3(H20). 



5.2. Choice of the Temperature Profiles 

Temperature profiles of dust cocoons around embedded massive stars remain one of 
the most difficult questions that need to be answered. Since the emission distribution is an 
interplay between the temperature and density distributions, the temperature profile directly 
affects the extraction of the density profile. The UV radiation from the protostar is absorbed 
by the dust cocoon and re-radiated as dust thermal emission, whose spectrum peaks in the 
JR. In general, the temperature profile refiects the diffusion of the IR photons. The innermost 
part of the core is very dense and optically thick to the IR photons while the outer part of 
the core should be optically thin. The transition between the two regimes possibly happens 
around r ~ 10^ AU. To date, observations have been limited by angular resolution and have 
mainly focused on the density profiles in the outer regime, where the temperature profile 
has a simple radial dependence of r~°'^. The temperature of the inner regime can only be 
calculated with radiative transfer models. However, simple theoretical calculations do not 
account for a wide variety of processes that may occur in accreting fiows around massive 
protostars, for example, radial variation of dust properties, effects of magnetic field, etc. 
Besides, the feedback of numerous chemical reactions to the heating and cooling processes 
has generally been neglected. Our dust continuum observations have partially resolved the 
inner region, where the temperature profile may differ significantly from that of the outer 
regime. We have chosen two extreme cases, optically thin and optically thick, to investigate 
the corresponding parameter spaces. 

In the outer regime, theoretical studies (Wolfire & Cassinelli 1986; Osorio, Lizano, 
& D'Alessio 1999) have shown that the temperature profile is expected to asymptotically 
approach a power-law index of g ^ 0.4 and is a weak function of the dust opacity law 
(Wolfire & Cassinelli 1986). Assuming that the dust thermal emission can be characterized 
by a modified Planck function of a single temperature To, we can write the local radiative 
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Fig. 3. — Continuum spectra of the two components in W3(H20). Dust thermal emission 
is detected at milUmeter wavelengths. Both components show rising spectra with spectral 
indices of 3.0 ± 0.4 and 2.9 ± 0.4, which are smaller than the spectral index of 3.6 from 
previous studies. At lower frequencies, the east component W3(H20)A shows synchrotron 
emission (Reid et al. 1995) while the west component W3(H20)C has a rising spectrum 
indicating another possible protostellar object (Wilner et al. 1999). 
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equilibrium equation as 

„ , , const 
K,B,{TD)du = (6) 

'0 ^ 

If the dust opacity law takes the form of K,y oc z/^, the integration gives a temperature 
distribution of To oc r^'^^^^^^^K Since surveys in star-forming regions usually measure the 
opacity law to be in the range of /3 = 1 — 2, the power- law index of the temperature profiles 
therefore falls into the range oi q — 0.4 — 0.33. The observed overall dust opacity law of 
Ki, oc v in W3(H20) will have q — 0.4 in the optically thin regime. 

In the inner regime where the optical depth is large for the IR photons, one can derive 
a simple relationship between the density and temperature distributions using the standard 
diffusion approximation (Osorio et al. 1999; Kcnyon, Calvet, & Hartmann 1993) 

L — — = const. (7) 

oKrP dr 

where L, p, ass are the luminosity, density, and Stcfan-Boltzmann constant, respectively. 
The Rosseland mean opacity kr{T) is an averaged opacity defined as 



1 as, , 

dv 



^ K„ dT 
— /•(- 



dB ■ 
-i^dv 
dT 

When the opacity law is k^, oc the Rosseland mean opacity can be easily calculated to be 
kr(T) oc T. Assuming a power-law dependence on radius for the density and temperature 
distributions, p oc r"^ and T oc r"'', one can apply the Rosseland mean opacity of kr{T) oc T 
to Eq. (7) and obtain a temperature profile with a power-law index oi q = {p + l)/3. Given 
the observed radial emission profile with a power-law index of approximately 2 = p -\- q 
(see Sect. 5.4), the largest possible value for q will be 0.75. Although the exact temperature 
behavior is more complicated than merely a single power-law, we have chosen the two extreme 
cases, q — 0.4 and 0.75, for our interpretation of the data. 



5.3. Emissivity Profiles 

Since the 1.4 mm beam (0'.'26) is about 5 times smaller than the extent of each source 
in W3(H20), continuum emission profiles can be calculated by taking averaged intensities 
in successive annuli centered at each continuum peak. Fig. 4 shows the emission profiles by 
plotting the averaged intensities against the outer radii of the annuli with a spacing of one 
synthesized beam to avoid undesired correlation between data points. In order to reduce 
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the interference between the two sources, we have masked one source while taking annular 
averages on the other source. The irregular extended emission to the northwest is also 
excluded from the emission profiles of source C. The uncertainty of each annular average, dj, 
is defined as the ratio of the root-mean-square (rms) in the map to the square root of the 
number of beams within the annulus of interest: 



where a is the rms of the map; Ni and A^Beam are the numbers of pixels in the ith annulus 
and one beam, respectively. 



The optimized solutions are hsted in Table 5 for the two selected q values and the results 
are shown in Fig. 4. Overall, the emission profiles have power-law indices of approximately 
-2 for both sources. The emissivity is proportional to the product of the temperature and 
density as e(r) oc pQTQr~''P^''\ If the temperature distribution follows T oc r""'^, the under- 
lying density profile will have a power-law index of p ^ 1.6. The density distributions may 
have an even smaller index of p ~ 1.2 if the temperature profile is as steep as T oc r~°'^^. 
In order to fit the innermost data point at Cy.'lT, we introduce a dust-free cavity of radius 
i?cav, which may correspond to the dust destruction radius. The value of it!cav is in the range 
of 200 to 250 AU, much smaller than the angular resolution of 0'.'26 (~ 570 AU). The real 
temperature distribution shall deviate greatly from a single power-law temperature profile, 
especially near the dust destruction radius. A reliable estimate of i?cav will need better 
resolution and more elaborate calculations. 

The self-similar collapse of an isothermal sphere (SIS) (Shu 1977) is able to explain the 
collapse of low-mass protostellar cores. This inside-out model describes the quiescent part 
of an isothermal cloud core with a density profile of p oc and propagates the collapse 
signal with an expansion wave, behind which a p oc r~^/^ law holds for the freely falling inner 
envelope. On the other hand, earlier observations of massive star-forming cores on larger 
scales suggested the mean density index p in a range between 1.6 to 1.8 with large scatters 
in the samples (van der Tak et al. 2000; Beuther et al. 2002; Mueller et al. 2002). The 
cores may be supported by turbulence and magnetic fields (Myers & Fuller 1992; McKee & 
Tan 2003). Studies of dense cores associated with massive stars, L1688 and L1204, (Myers 
& Fuller 1992) suggested that the support from turbulence or magnetic fields most likely 
dissipates on the scale of 0.01 pc, where we observed the density profiles of the two sources 



a 



(9) 




5.4. Density Profiles 
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Fig. 4. — The 1.4 mm continuum emission profiles of annular averaged intensities versus 
the outer radii of the annuli. The annular averages are taken with a spacing of the full 
synthesized beam so the data points are largely independent. The optimized solution for 
q — 0.4 (solid curve) and that for q — 0.75 (dashed curve) fit the data equally well. The 
cutoffs at large radii are caused by the finite sizes of the dust distributions. The turnovers at 
the innermost radii result from a combined effect of the central cavities and the convolution 
with the beam. 
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Table 5. Optimized Fits for the Continuum Models 



Parameter 




0.4 




0.75 


Source 


A 


C 


A 


C 


A mm 


1.00 


1.14 


1.00 


1.15 


V 


1.52 ± 0.17 


1.52 ± 0.27 


1.15 ± 0.16 


1.14 ± 0.27 


po (10-18 gcm-3) 


3.59 ± 0.42 


2.78 ± 0.47 


3.48 ± 0.36 


2.48 ± 0.41 


Rd (10=^ AU) 


2.05 ± 0.09 


2.06 ± 0.16 


2.05 ± 0.09 


2.06 ± 0.16 


i?eav (10' AU) 


2.48 ± 0.51 


2.17 ± 0.83 


2.50 ± 0.46 


2.18 ± 0.81 


Parameters set by the CH3CN models 






^0 (K) 


197 


172 


205 


192 


Parameters derived from models: 










0.094 


0.075 


0.092 


0.071 


Ma (Mo) 


0.043 


0.034 


0.057 


0.041 


M, (Mo) b 


4.3 


3.4 


5.7 


4.1 


(10^ cm-3) be 


1.07 


0.83 


1.04 


0.74 



^Averaged optical depth, (r), within O'.'S around the peak. 

'^A gas-to-dust ratio of 100 is assumed. 

'^The number densities of H2 at Tq, calculated from po- 
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in W3(H20). The small spatial extent of our density profiles may favor the explanation of 
free-fall collapse. Therefore, we suggest that both sources in W3(H20) are in an active phase 
of accreting material from their surroundings. 

In order to estimate the optical depth of each clump, we calculate an averaged optical 
depth, (r), which is weighted by the solid angle of each annulus, 

1 - e<^> = ^ — jr— ^- , 10 

J27rada ^ ' 

where a is the angular distance from a continuum peak. The integral is performed within 
0'.'5 around each peak, the same region where we integrate the continuum flux density. In 
general, the averaged optical depths are in the range of 0.07 to 0.09. Considering a clump of 
uniform density and temperature with (r) = 0.09, the error of the 1.4 mm flux density made 
by the optically thin approximation is (t)/(1 — e"^"^^) ~ 5%. This in turn contributes an 
additional error of only 0.07 in the spectral index. Therefore, the optically thin assumption 
should be reasonable and the small spectral indices are not the result of large optical depths 
at 1.4 mm. 



5.5. Mass Estimates 

Table 5 also lists the estimated masses and maximum optical depths of the dusty nebula 
around each source. The mass of the dust, M^, can be computed simply by integrating the 
density over the volume of the clump. Once the mass of the dust is known, the clump (gas) 
mass, Mg, can be obtained by assuming a conventional gas-to-dust ratio of MgjMd = 100. 
Over all, the nebular mass of source A is in the range of 4 to 6 Mq while that of source C is in 
the range of 3 to 4 Mq. Nevertheless, these masses exclude any mass contained in the central 
compact objects. The number densities of H2, nu^, at the reference radius, tq = 500 AU, 
are also listed in Table 5. In general, the number densities are roughly 10^ cm~'^, which are 
much higher than riu^ — 2 x 10^ cm~^ as suggested by Helmich & van Dishoeck (1997) in a 
study of single-dish spectral line surveys. 

A gas-to-dust ratio of Mg/Md ~ 100 (Savage & Mathis 1979) is the most commonly 
accepted value, where it was derived from comparing the visual extinction with the hydrogen 
column density, and should be good within a factor of 2 (Hildebrand 1983). This ratio 
serves as a coefficient to convert dust mass to gas (clump) mass and is used in the mass 
estimates of Table 3. The largest uncertainty in the density (or mass) determination comes 
from the chosen dust opacity, where various estimates of fi:i.4mm span a range from 0.2 to 
3cm^g~^ (Henning, Michel, & Stognienko 1995; Beckwith & Sargent 1991). Observational 
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determinations of mass opacity require a few assumptions, each of which brings uncertainties 
into the final result (Hildebrand 1983). In addition, various environments may carry dust 
grains of different types, which also contribute to the variation in the derived mass opacity. 
On the other hand, theoretical calculations suffer from unknown dust shape and chemical 
composition. A mixture of small graphite and silicon spheroids may have a mass opacity as 
low as fi;i.4mm — 0.2 cm^g~^ (Draine & Lee 1984) while grain growth with fractal dimension 
can increase Ki.4mm by a factor up to 30 (Wright 1987). Based on a more sophisticated model 
of the coagulation process (Ossenkopf 1993), Ossenkopf & Henning (1994) systematically 
computed and tabulated dust opacities in dense protostellar cores. They suggested a dust 
opacity of 1.11 cm^ g~^ for a core of riu^ — 10^ cm~^ with non-negligible warming by near-IR 
radiation. 

The opacity of Ki.4mm = 1-84 cm^ g~^ in our continuum models is calculated by applying 
the observed power-law indices of /3 ~ 1.0 while maintaining the normalization at 250 /im: 



where the value of 10 cm^g~^ at 250 /im is given by Hildebrand (1983). This is a similar 
approach to what has been suggested by Beckwith & Sargent (1991) for deriving the masses 
of circumstellar disks around low-mass young stars. The dust opacity of 1.84 cm^g^^ is a 
factor of 1.66 higher than the computed value of 1.11 cm^g^^ (Ossenkopf & Henning 1994) 
and may underestimate the nebular masses by 40% if one adopts the computed opacity. 



A temperature measurement is important to estimate the luminosity of an embedded 
protostellar object as well as to break the degeneracy of po^o in our continuum models. 
Temperature can be derived when multiple transitions of a molecular species with a great 
range of excitations are observed. Symmetric-top molecules with large dipole moments, 
such as methyl cyanide (CH3CN), are ideal "thermometers" to probe the gas temperature. 
This type of molecule has groups of hues with greatly different excitation energies but close 
in frequency, so that the spectral correlators can be adjusted to observe these transitions 
simultaneously. This reduces the uncertainty from calibration, antenna gains, Mv-sampling, 
etc. In addition, the radial velocity can be determined with high precision by measuring 
the Doppler shifts in all the transitions. Line observations with good spectral resolution are 
necessary to obtain kinematical clues for understanding the relationship of the two continuum 
sources. 
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6. METHYL CYANIDE LINE OBSERVATION RESULTS 
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6.1. Kinematics of the bineiry 

Assuming that all the K components have the same velocity, we can simultaneously fit 
Gaussians to each of the five K components to trace the bulk motion of the source. The 
results of our Gaussian fits are shown as the green curves in Fig. 5. Only data above the 
2(7 level (yellow curves) are included in the fits. The velocities found in Gaussian fits are 
Vr — —51.37 ± 0.05 kms~^ at source A and Vr — —48.56 ± 0.05 kms~^ at source C. There- 
fore, the velocity difference between the two continuum peaks is l^Vobs — 2.81 ± 0.10 kms~^. 

Consider binary motion consisting of two circular orbits at an inclination i. The separa- 
tion, a, of the binary members is the sum of the semi-major axes of the two orbits. Let the 
relative velocity between the binary be — Vi — V2. At the maximum projected separation, 
a, the radial velocity difference is also the largest, l^Vr — sin i. The total mass of the 
binary can be expressed as 

Mbina^y - Mi + = Z,. , (12) 

G sm I 

where G — 6.673 x 10~^ g~^ cm^ s~^ is the gravitational constant. Suppose that the two 
sources in the hot core happened to be at their maximum projected separation. The binary 
separation is just the observed separation, a — aobs — 2.43 x 10^ AU, and the radial velocity 
difference is — ^Vobs — 2.81 ± 0.10 kms~^. Since the trigonometric function in Eq. (12) 
is always less than unity, the minimum binary mass will be 

Mw3(H.o) > - 1.13 (^^-^ j [-WPJj) ^® 

= 22 M0. (13) 

This sets a lower limit of 22 M0 for the total mass of the two protostars. The radial velocity 
difference that we observed here is in a reasonable range for binary motion. Therefore, we 
suggest that the two continuum peaks are members of a binary and the velocity difference 
is due to the binaries orbiting around each other. 

Sutton et al. (2004) have also observed a small radial velocity difference of 1.6 kms^^ 
between source A and C with several methanol lines at low spatial resolution. They suggested 
that the velocity difference is caused by binary motion. The smaller velocity difference should 
be a result of averaging over a larger solid angle. Besides, methanol emission seems to emerge 
from a region larger than that of methyl cyanide. The geometry and kinematical structures 
of the circumstellar envelopes will also affect the observed radial velocities. 
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Fig. 5. — Spectra of methyl cyanide emission observed at source A and source B (shown in 
histograms). Only data points above the 2a level (yellow curves) are included for further 
analysis. The radial velocities, t;^, are found by simultaneously fitting 5 Gaussians (green 
curves) ; each Gaussian has its own width and amplitude but the five K components have the 
same radial velocity. The optimized solutions for q = 0.4 are shown in blue dashed curves 
while those for g = 0.75 are shown in magenta dashed curves. The big K values indicate the 
K components of the CH3CN J = 12 — > 11 transitions and the small k values indicate the 
corresponding K components of its isotopomer CH3 ^^CN. 
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6.2. Boltzmann Population of Levels 



The methyl cyanide spectra (Fig. 5) clearly show a suppression in the K = 3 component. 
The ratio of the peak antenna temperature of A' = 2 to that of i^T = 3 is rougly unity, 
which is much smaller than the expected optically thin ratio of roughly 2 based on the 
product of the statistical weight and line strength, QirSik (Table 3). It is evident that 
the K — 2) component is optically thick and its antenna temperature should be close to 
the kinetic temperature of the gas if the emission fills the beam. On the other hand, the 
K — Q component with an upper energy level of E'up = 326 K is also detected at both 
sources, suggesting the presence of hot gas components. Nonetheless, the observed antenna 
temperature of the K = ?> line is only T4 ~ 60 K, which is insufficient to explain the detection 
of the K = Q component. Therefore, the methyl cyanide emission must come from a region 
significantly smaller than our 1" beam, equivalent to a radius of 10^ AU. 

Although the emission is spatially unresolved, the antenna temperature of the optically 
thick lines can still be used as a lower limit of the kinetic temperature of H2- Given this 
minimum value of 60 K, we can show that the observed methyl cyanide lines satisfy the 
conditions for LTE, which will give the Boltzmann level population without further assump- 
tions. The coUision rate coefficients, Kui, of the observed lines are only involved with the 
fimdamental rates Q{L,M = 0) (see Sect. A). In the case of CH3CN (12, X) (ll,i^) 
transitions, K^i can be estimated by its largest leading term in Eq. (A3) 



where the 3 — j symbol is about 1.5 for all the K components. At a temperature of 
60 K, Q(1,0) = 3.22 x 10"^° cm^s'^ so the collision rate coefficients will be K^i > 1.67 x 
10~* cm^s"^. Helmich & van Dishoeck (1997) estimated a molecular hydrogen density of 
riHa = 2 X 10^ cm~'^ towards W3(H20). Therefore, the minimum value for downward collision 
rate will be Cui ^ 3.4 x 10~^ s~^. This minimum rate of Cui will be even faster if the number 
densities of — 10® cm~^ in Table 5 are used. On the other hand, the maximum value 
of the spontaneous emission rates among the observed fines is = 9 x 10~^ s"^ from the 
K — 2 component (Table 3). In other words, the collisional de-excitation rates are faster 
than the spontaneous emission rates by at least a factor of 40. The transitions that we 
consider here should have level populations that follow the Boltzmann distribution at the 
kinetic temperature of H2. 
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6.3. Model Results 

The velocity of each source has been obtained in the kinematics study (Sect. 6.1). We 
have assumed that methyl cyanide density has the same power-law dependence on radius 
as that of the dust. However, the spatial extent of the methyl cyanide emission, -RchsCN; is 
smaller than that of the continuum emission, i?cont- Since the extent of the methyl cyanide 
is unknown, the emission size, i^cHsCN, is a parameter determined by the chi-square (x^) 
minimization. The optimized solutions for the two selected q values are listed in Table 6 
and the results are shown in Fig. 5. The values for source A are smaller than those for 
source C. In general, the values of q do not dramatically change the goodness of the fits. 

Over all, the Gaussian line width h is much larger than what one expects in the case 
of merely thermal motion, hth — {kT/mcnaCNY^^ — 0.07 kms~^ at T = 200 K. A velocity 
dispersion greater than the thermal velocity dispersion has been observed in low-mass dense 
cores within the scale at which turbulence dissipation takes place (Goodman et al. 1998). 
However, the thermal component of the line width comprises only an insignificant fraction 
of the total line width in our case. It is not clear what mechanism is responsible for the 
nonthermal component: outfiow, rotation, infall collapse, or small-scale turbulence may all 
contribute to the line width. Besides, the velocity gradient between source A and G or the 
overlap of the two clumps may also increase the line widths. 

The line width at source G (3.5 kms~^) is much larger than that at source A (2.8 kms~^), 
suggesting an active, embedded source. Ghemical differentiation has also been observed: 
Wyrowski et al. (1999) have detected the i>7 = 2, J = 24 — 23 transition of HG3N (excitation 
at 722 K above the ground) only at source G and the G2H5GN (242,22 — 232,2i) emission 
peaked at source G. Our spectra also show a stronger G2H5GN (252,24 — 242,23) emission at 
source G. This additional evidence also supports source G being a more active source rather 
than just an over-dense clump. 

The GH3 ^^GN lines can be used to constrain the optical depths of the main isotopomer 
lines. Gomparing the measured antenna temperatures of the K — 2 components of the two 
isotopomcrs, the GH3GN line is brighter than the GH3 ^^GN line by roughly a factor of 4.4. 
The intensity ratio between the two species can be approximated by (CH3 ^^CN/CHaCN) = 
{^^C/^'^C) [t/ (1 — e~'^)], where r is the optical depth of the CH3CN line. Assuming a distance 
of Dqc = 10 kpc to the Galactic Center, the carbon isotope ratio is about (^^GO/^^GO) = 
82.6 in the vicinity of W3(0H) (Wilson & Rood 1994). Therefore, the optical depth of the 
K — 2 component of the GH3GN fine should be less than r < 19. 
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6.4. Mass Estimates 

The virial mass, Mvu-, is a function of the hne width, b, the power-law index, p, of the 
density distribution, and the size of methyl cyanide emission, i^cHsCN (see Appendix B): 

On average, the virial mass of source A is about 12 Mq while the mass of source C is 
about 28 Mq (Table 6). The derived virial masses are large due to the large line widths 
and should be considered as upper limits for the mass. A massive star-forming core is an 
entity of complicated kinematics so a dynamical equilibrium most likely has not yet been 
achieved considering the youth of a hot core. The virial theorem is based on two fundamental 
assumptions: first, the objects form a relaxed system whose motion is gravitationally bound; 
second, the lines are not significantly broadened by other effects, such as large optical depth, 
systematic velocity gradients, etc. Any deviation from the above two conditions will increase 
the value of the derived virial mass. 



6.5. Luminosity Estimates 

The simplest approach for calculating the luminosity is to model the emission as a black 
body like that from a stellar surface: 

L = ATrr^asBT^ 

where asB = 5.67 x 10~^ ergcm~^s~^ is the Stefan-Boltzmann constant. This will give 
only a rough estimate of the luminosity. The luminosity for source A is in the range of 
(1.5 — 1.8) X 10^ Lq while that of source C is in the range of (0.9 — 1.4) x 10^ Lq, equivalent 
to zero-age-main sequence (ZAMS) spectral type BO. 5 - BO stars (Panagia 1973). Previous 
studies have also suggested a total luminosity of 10'^ Lq for the entire hot core (Turner & 
Welch 1984; Wilner et al. 1995; Wyrowski et al. 1999). The stellar masses of ZAMS spectral 
type BO. 5 - BO stars are in the range of 11 to 15 Mq (Hanson, Howarth, & Conti 1997). 
Therefore, the total mass of the binary system is in the range of 22 to 30 Mq, in good 
agreement with the minimum kinematical mass of 22 Mq (Sec. 6.1). 
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An interesting question that arises with high luminosity is whether the radiation pressure 
acting on dust grains is large enough to halt the infall (Larson & Starrfield 1971; Kahn 1974). 
Within factors of order unity, radiation pressure can reverse spherical infall if the luminosity 
to mass ratio, L/M, of the protostar exceeds a critical value of 700 Lq/M©. This critical 
value, however, varies with the assumed dust properties and the geometry of infall processes 
(see, e.g., the review by Shu, Adams, & Lizano (1987)). Over all, the ratios of luminosity 
to nebular mass, L/Mg, in our model results (Table 6) are larger than this critical value but 
within a factor of 5. The masses contained in the central protostars also add uncertainties 
to the observed ratios. 

When the expansion of an H II region is considered, the opacity is mainly due to the 
scattering of free electrons, and the Eddington limit gives a luminosity to mass ratio of 
L/M — 3.2 X 10^ L0/M0. Our luminosity to nebular mass ratios, L/Mg, are smaller than 
the Eddington limit by roughly one order of magnitude. 

6.6. Thermal Coupling between the Dust and Gas Components 

The dust warmed by radiation from embedded sources determines the thermal structure 
of hot cores. Since line heating by atoms and molecules is significantly weaker than the 
continuum absorption by the dust, the gas absorbs very little of the stellar radiation and is 
further shielded by the dust. Therefore, the dust is hotter than the gas and is the dominant 
heating agent for the gas. Two possible heating processes transfer energy to the gas: collisions 
with dust grains and line absorption of the IR radiation field. Our goal here is to show 
whether thermal coupling between the gas and dust components can be achieved in hot 
cores rather than to make detailed calculations. 

The thermal emission of the warm dust builds an intensive IR radiation field character- 
ized by a radiation temperature equal to the dust temperature, T^. The molecular hydrogen 
has a great number of vibration-rotation transitions {e.g. fundamental f = ^ 1 lines at 
2.4 yum) and pure rotational lines {e.g. J = 2 ^ at 28 /xm) in the IR band that can respond 
to the ambient IR radiation field. If the optical depths of these transitions are large in the hot 
core, the hydrogen excitation temperature should approach the radiation temperature at 
all time. Then collisions between molecules will convert the internal energy to (macroscopic) 
kinetic energy. On the other hand, if the optical depths are small, line heating will be less 
important for warming up the gas. We have estimated the optical depths for the transitions 
of {vJ — > v' .J') = (11 — > 01) at 2 fj,m and i' = 0J = 2— >0at28 fim. Given a minimum 
line width of 2.7 kms~^ and a maximum column density of A^H2 — 4.4 x 10^^ cm~^ in our 
models, the maximum optical depth of H2 transitions in the bands of 2 /^m and 28 /xm is 
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Table 6. Optimized Fits for the Methyl Cyanide Line Models 



Parameter 




0.4 




0.75 


Source 


A 


C 


A 


C 


■ 

A mm 


1.84 


2.46 


1.86 


2.50 


no (cm-^) 


2.70 ± 0.28 


2.00 ± 0.14 


2.74 ± 0.25 


1.95 ± 0.14 


To (10^ K) 


1.97 ± 0.14 


1.72 ± 0.09 


2.05 ± 0.13 


1.92 ± 0.08 


b (MHz) 


1.98 ± 0.05 


2.68 ± 0.06 


2.00 ± 0.05 


2.70 ± 0.06 


i^CHsCN (10^ AU) 


7.16 ± 0.27 


9.01 ± 0.30 


7.13 ± 0.18 


8.88 ± 0.24 


Parameters set by the continuum models: 






P 


1.52 


1.52 


1.15 


1.14 


i?eav (10' AU) 


2.48 


2.17 


2.50 


2.18 


Parameters derived from the models: 






b (kms"^) ^ 


2.69 ± 0.07 


3.64 ± 0.08 


2.72 ± 0.07 


3.67 ± 0.08 


Mvir ( Mq) 


11.7 


26.9 


13.1 


29.7 


L (10^ Lo) 


1.5 


0.9 


1.8 


1.4 


L/Mg (10^ Lq/Mq) 


3.5 


2.6 


3.2 


3.4 



^Velocity line widths are calculated from the frequency line widths. 
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about 1. The optical depths of H2 transitions are too small to ensure good thermal couphng 
through line heating. 

Dust-grain collisions are believed to be an important heating process to explain the 
gas temperature of the molecular clouds with embedded sources (Goldreich & Kwan 1974; 
Takahashi, Silk, & HoUenbach 1983). The heating rate per H2 molecule by collisions with 
dust grains is given by Takahashi et al. (1983) as 

= nu^dva2k{Ta-Tg) (18) 



- i.4xio-„„W-r4iL)(^) 

ergs-^Hs"' (19) 



1.0 X 10-21 cm2 

where nn is the number density of hydrogen nuclei, T^d — ndcrd/n-H. is the projected grain 
area per H nucleus, and a is the accommodation coefficient of the grain surface for an H2 
molecule (Burke & HoUenbach 1983). The mean speed of H2, v = {SkTg/TTmn^Y^'^ , is about 
1 kms"^ at gas temperature Tg ~ 100 K. Although the mean speed increases with the gas 
temperature, the value should be good for an order-of-magnitude estimate. The gas heating 
rate is roughly equal to rgrain/'^H2 = c^lf^l^s ~ Td)]/dt. Therefore, the gas-dust relaxation 
time, tgd, is approximately 

2-5 X 10^^ 

tgd = s. (20) 

Given the minimum dust density of po — 2.48 x 10"^^ gcm"^ (see Table 5) in our models, 
the corresponding number density of H nuclei is nn = 1.5 x 10^ cm"^ if a gas-to-dust ratio of 
100 is assumed. Therefore, the gas-dust relaxation time is about tgd ^ 1.7 x 10^ s = 5.4 yr, 
which is significantly shorter than the lifetime of hot cores. As long as the stellar radiation 
does not vary greatly on the time scale of the gas-dust relaxation time, the gas and dust 
components are well coupled thermally by collisions. 



6.7. The Nature of the Velocity Gradient 

The linear velocity gradient suggested by Wyrowski et al. (1997) is also detected in our 
methyl cyanide lines. For the purpose of easy comparison, we only perform the Gaussian fits 
on the K = 3 component and the result is shown in the left panel of Fig. 6. The position 
uncertainty, ae, of a Gaussian emission feature due to thermal noise is given by Reid et al. 
(1988) as 

/ 4V''^^FWHM 1 „ ^FWHM 

"'^W VlEl SNR = "-"^ W 
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where ^fwhm is the full-width at half-maximum of the emission feature, and SNR is the 
signal-to-noise ratio, which is set equal to the peak intensity divided by the rms noise well 
away from the emission. 

A similar velocity gradient along the East- West direction across W3(H20) is clearly 
seen. With improved angular resolution, our data are able to trace over a broader velocity 
range with better sampling and provide a few important details. Our position drift shows 
turnovers, where the peak positions stay unchanged at the two ends, and a twitch, where the 
velocity has a larger increment (from —52 kms^-*^ to —50.5 kms~"^) at 5'.'5 east to the phase 
center. These irregular behaviors in the velocity gradient can result from two unresolved 
clumps with a small radial velocity difference. In order to inspect this hypothesis, we have 
generated a binary data cube by adding linearly the methyl cyanide models of the two 
sources and performed the same Gaussian fitting to the model images. The model images 
show similar behaviors (right panel of Fig. 6) and suggest that the binary rotation can explain 
the irregularities in the velocity gradient. To have a better understanding, we have taken a 
velocity-position slice that passes through the two continuum peaks (Fig. 7). At first sight, a 
velocity gradient seems obvious; however, the sharp edges suggest that the velocity gradient 
is a combination of two unresolved clumps. The same position-velocity slice with the model 
data cube shows that two connected ellipses are expected from a partially resolved binary 
motion. 

Our angular resolution of 1" has improved by nearly a factor of 2 compared with the 
previous study (Wyrowski ct al. 1997). This better angular resolution has provided impor- 
tant clues to reveal the rotation of an unresolved binary. When the emission is only resolved 
in one direction, a Gaussian function can be a misleading approximation for the emission 
morphology. We conclude that the linear velocity gradient is actually a result of unresolved 
binary rotation rather than the Keplerian rotation of a disk or the expansion of a bipolar 
outflow. 



7. SUMMARY 

1. The W3(H20) hot core is the best case to date that shows a massive protobinary 
system at such an early evolutionary stage that it is still associated with nebular 
emission from the surrounding dust cocoon. This hot core is also a good example 
showing fragmentation as part of the star formation process for massive stars. 

2. We have obtained the best angular resolution that has ever been observed for W3(H20) 
at millimeter wavelengths: the synthesized beam size reaches 0'/26 at 1.4 mm and 0'/4 
at 2.8 mm. Our continuum maps show a double-peaked morphology, suggesting an 
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Fig. 6. — Right ascension offset of the Gaussian peak position versus the radial velocity in 
W3(H20). A velocity gradient is clearly seen in the observed K — 3 component (left) and 
the model images (right) . The velocity gradient generally agrees with the result in Wyrowski 
et al. (1997). The discontinuity of the position drift on the edges and the twitch at offset 
5'/4 suggest the presence of multiple clumps. 
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POSITION POSITION 

Fig. 7. — Velocity-position plot of the K = 3 component of the CH3CN J = 12 — 11 
transition (left) and the corresponding plot for the model data cube (right). Contour levels 
correspond to 5.7 to 51.3 by 5.7 {2a) K. 
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embedded protobinary system with a projected separation of 2.43 x 10^ AU. 

3. The kinematics gives a radial velocity difference of 2.8 kms~^, corresponding to a 
minimum binary mass of 22 Mq. The continuum models provide a nebular mass of 4 
to 6 Mq for source A and 3 to 4 Mq for source C. 

4. The two binary members have a dust opacity law of k^, oc u, suggesting that grain 
growth has begun. 

5. If the temperature profile follows T oc r~'^-^, the density distribution for each source has 
a power-law index of —1.52, which is close to the power-law index of —1.5 for free-fall 
collapse. 

6. The luminosity of source A is in the range of (1.5 — 1.8) x 10^ Lq and that of source C 
is in the range of (0.9 — 1.4) x 10^ Lq. The luminosities of the two binary members 
correspond to ZAMS spectral type BO. 5 - BO. Since a star of spectral type earlier than 
B3 will produce an observable Hll region, we expect each member of the binary to 
develop its own H ll region in the future. 

7. Based on the brightness temperature argument of the five K components, the methyl 
cyanide emission should emerge from regions significantly smaller than the extent of 
the warm dust. The methyl cyanide line model provides a temperature estimate of 
200 K for source A and 182 K for source C at the reference radius of 500 AU. 

8. Excluding any mass contained in central protostars, the derived luminosity to nebular 
mass ratio, L/Mg, is larger than 700 Lq/Mq, but smaller than the Eddington limit by 
one order of magnitude. 

9. The virial masses derived from the line widths set upper limits for the masses at 12 Mq 
for source A and 28 Mq for source C. 

We wish to thank Mel Wright for help with the BIMA data analysis and plotting pro- 
grams, MIRIAD and WIP. Thanks to Dick Plambeck for help with many questions that 
arose during this research. Thanks also due to Al Glassgold for the discussion on the ther- 
mal coupling between the dust and the gas. This work was supported in part by NSF grant 
AST93- 14847. 



A. PROPERTIES OF METHYL CYANIDE MOLECULES 

The energy levels of a symmetric-top molecule can be described by two quantum num- 
bers, J and K, where J is the total angular momentum and K is its projection on the 
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symmetry axis of the molecule. The rotational energy Ejk of a symmetric-top rotor includ- 
ing centrifugal stretching can be written as (Sutton et al. 1986; Anttila et al. 1993) 



where A, B, Dj, Djk, and Dk are molecular constants and we have adopted the values 
given by Pearson & Mueller (1996). The energy level diagram of methyl cyanide is showed 
in Fig. 8. Since a symmetric-top molecule has no dipole moment perpendicular to the 
symmetry axis, a radiative transition can not change the angular momentum along this axis. 
Therefore, radiative transitions are only allowed between adjacent J levels within a single K 
ladder. The spacing between adjacent J levels only decreases slightly in different K ladders. 
As a result, the transitions of adjacent J levels, e.g. J ^ J — 1, in various K ladders 
happen at nearly the same frequency and can be observed simultaneously. Each successive 
K component is shifted to a slightly lower frequency due to the centrifugal distortion. 

The three identical hydrogen nuclei in methyl cyanide create threefold symmetry. When 
a pair of hydrogen nuclei are exchanged, the combined effect of 120° rotation and nuclear 
spins brings additional symmetry requirements that further divide the molecules into two 
distinct species: ortho- and para-CHsCN. Rotational levels with K = 3n {n is an integer), 
belong to the ortho-species, while levels with K = 3n ± 1 belong to the para-species. The 
ortho- and para-species are sometimes referred to as the A and E species, respectively. The A 
and E species are essentially independent since the probability of changing nuclear spin state 
in a collision is very small and radiative transition between species is forbidden. Because of 
the symmetry requirement, the statistical weight, qjk, of the A species is twice larger than 
that of the E species with the exception of the K — ladder, where half of the levels are 
forbidden due to the total cancellation of the wave function. 

The coUision rate coefficients for downward transitions, K^i, of methyl cyanide were 
first calculated by Cummins ct al. (1983) and expanded to higher energy levels and higher 
temperature up to 140 K by Green (1986) to better match the conditions in star-forming 
regions. The basic idea is to express collision rate coefficients in terms of a few fundamental 
rates, Q{L, M), that are simply related to excitation out of the lowest J = K = level. The 
collision rate coefficient for the transition (J, K) — > (J', K') of a symmetric-top molecular 
can be expressed as 



Ejk 
IT 



BJ{J+1) + {A-B)K^ 

-Dj J\J + If - DjK J(J + 1)K^ - Dk 
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Fig. 8. — Energy level diagram of CH3CN for J levels up to J = 12 and K ladders from 
K — through K — 6. The observed transitions are indicated by arrows. The A and E 
letters label the two independent species of K ladders distinguished by symmetry. 
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where the large bracket is a 3-j symbol. The second term in the large square bracket on the 
right-hand side of Eq. (A3) is much smaller than the first term and was neglected in our 
calculations in Sect. 6.2. The fundamental rates, Q{L, M), that we used in this study are as 
given by Green (1986). 



B. VIRIAL MASS FOR A POWER-LAW DENSITY DISTRIBUTION 

Given a power-law density distribution of p = Po{r / rc,)^^ . the enclosed mass within a 
radius of r is m(r) = An po (3 —p) . The gravitational potential energy, U, for a clump 

of radius R can be calculated by 

rR 

U — —AnG i ra{r) p{r)rdr 



JO 

(3-p)(5-2p) 
= 1^^, (Bl) 

(B2) 

where M — m{R) = Att po / (3 — p) is the total mass of the clump. The kinetic energy 

can be written as = 3NkT/2 = 3Mb'^/A, where is the total number of particles, and 
b is the gaussian linewidth in Eq. (5). Assuming the gas is in quasistatic equilibrium, the 
gravitational potential energy of the clump can be related to its kinetic energy by the virial 
theorem, U = —2K. Therefore, the virial mass can be written as 
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